Interference Relay Channels - Part I: 
Transmission Rates 

Brice Djeumou, Student Member, IEEE, Elena Veronica Belmega, Student Member, IEEE, 
§ ■ and Samson Lasaulce, Member, IEEE 

(N : 
<\ 

^ ' Abstract 

, We analyze the performance of a system composed of two interfering point-to-point links where the transmitters 

^ ■ can exploit a common relay to improve their individual transmission rate. When the relay uses the amplify-and-forward 

protocol we prove that it is not always optimal (in some sense defined later on) to exploit all the relay transmit power 

"sj i and derive the corresponding optimal amplification factor. For the case of the decode-and-forward protocol, already 
^ '. 

, investigated in HI, we show that this protocol, through the cooperation degree between each transmitter and the relay, 
00 . 

lO ' is the only one that naturally introduces a game between the transmitters. For the estimate-and-forward protocol, we 

(N ■ 

' ' derive two rate regions for the general case of discrete interference relay channels (IRCs) and speciauze these results to 

^ ! 

' obtain the Gaussian case; these regions correspond to two compression schemes at the relay, having different resolution 

, levels. These schemes are compared analytically in some special cases. All the results mentioned are illustrated by 

■ simulations, given in this part, and exploited to study power allocation games in multi-band IRCs in the second part 

' of this two-part paper. 

Index Terms 

Amplify-and-forward, decode-and-forward, estimate-and-forward, interference relay channel, relay channel, re- 
laying protocols. Shannon transmission rates. 



I. Introduction 

The case of unlicensed frequency bands has become more and more important over the last decade in part 
because of the resounding success of Wi-Fi systems. In these bands, wireless devices generally communicate 
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in a non-coordinated manner which therefore leads to a scenario where signals interfere. In order to limit 
interference, devices have to transmit with a relatively (e.g., in comparison with a cellular phone) low power 
(17 — 20 dBm typically) which implies, in particular, some limitations on the communication range and rate. 
A quite natural way of improving the range, transmission rate or/and quality of the communications is to 
add relaying nodes that can be exploited by certain/all devices operating in the same frequency band. The 
mentioned scenario therefore gives one (but not the unique) strong motivation for studying the following 
system: a network comprising two transmitters communicating with their respective receivers and a relaying 
node that can be used by both transmitters. In this two-part paper, we analyze two key aspects of this type 
of networks, which is modeled by an interference relay channel (IRC) as defined in |[il|l2l and re-defined 
later on. In the first part, we focus on the relaying protocols for such a system. Our goal is to know how 
the three dominant classes of protocols, namely amplify-and-forward (AF), decode-and-forward (DF) and 
estimate-and-forward (EF), can be adapted, used and compared in terms of transmission rate in the context 
of interest. In the second part, we assume that the two transmitter-receiver pairs can communicate over 
several bands for which a relay is available on each of them (parallel IRC). The considered framework 
corresponds to the powerful paradigm of cognitive wireless networks Il3lll4l where transmitters can sense 
their environment and react accordingly. More specifically, in our context, the transmitters have to decide 
by themselves the best power allocation (PA) policy between the available bands in order to maximize their 
transmission rates. In order to analyze how the transmitters compete for the additional spectral resources and 
relaying nodes, we will exploit the results derived in the first part of the paper and also game theoretical tools 
which will be used to model the interaction between the transmitters and prove the existence of predictable 
and stable states (in the sense of Nash |[5l) for these networks. 

In this part of the paper, one of our goals is to assess the performance of the network under study when 
terminals implement good channel codes. This is why we will consider Shannon transmission rates as the 
performance measure. Note that, even in real systems where such good codes are not always implemented, 
considering Shannon rates can bring some insights on how to design the system. Indeed, there is a direct 
relationship between the achievable transmission rate of a user and its signal-to-interference plus noise ratio 
(SINR). Therefore, optimizing Shannon rates can also allow one to optimize performance metrics like the 
SINR or related quantities of the same type (e.g., the carrier-to-interference ratio). As far as the interference 
relay channel is concerned, note that the general idea of introducing a relay in a multiuser channel is 
quite natural and therefore has been exploited by other authors for other types of networks. For example. 
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the authors of H introduced the multiple access relay channel for which two transmitters can exploit a 
relay to communicate with a common receiver. The authors of |I3|[8l studied the dual of this problem, the 
broadcast relay channel where a single transmitter can cooperate with a relay to communicate with two 
receivers. As for the IRC itself, it has been introduced by [[Il[l2l- In [HI and O, the authors assumed that 
the relay implements the DF protocol and a (non-causal) dirty paper coding (DPC) scheme respectively, and 
derived achievable transmission rate regions for Gaussian IRCs. Part I of the present paper precisely aims at 
providing new results for the channel introduced by QO. As mentioned in the abstract, our contribution is 
essentially fourfold: 1) we derive an achievable transmission rate region for Gaussian IRCs when the relay 
implements a zero-delay scalar AF (ZDSAF) protocol and prove that there is an optimum amplification 
factor which does not always correspond to saturating the transmit power constraint at the relay (Sec. Hill ): 
2) assuming the DF protocol, we prove the existence of a game in Gaussian IRCs where the strategy of a 
transmitter is merely the correlation degree between its signal and the cooperative signal sent by the relay 
(Sec. HvT) : 3) we derive two achievable rate regions in the general case of discrete IRCs, corresponding to 
two EF-based protocols (Sec. IV]) using respectively one and two resolution levels (the Gaussian case readily 
follows); 4) we further discuss and compare all of these protocols by simulations (Sec. IVTl) . Before tackling 
these issues, we first present, in the following section (Sec. HH), the system model under consideration and 
main corresponding notations. 

II. System model and main notations 

Fig. [H depicts the system model considered in the whole paper except for the case of discrete input discrete 
output IRCs which will be studied in Sec. IV-AI The system comprises two transmitters, two receivers and 
a relaying node. The transmitters or information source nodes are denoted by Si and ^2, the relay node 
by TZ and the two receivers or destination nodes by Vi and V2. As indicated in Fig. [H the corresponding 
transmitted and received signals are respectively denoted by Xi, X2, X^, Yi, Y2 and Y^. The signal Xi 
(resp. X2) conveys the private message Wi (resp. W2) that Si (resp. ^2) sends to Vi (resp. I?2)- The general 
form of the received baseband signals is as follows: 

Yi = hiiXi + /i2i^2 + hj-iX,. + Zi 

< ¥2 = h22X2 + huXi + hr2Xr + Z2 (1) 
Yj. = hirXi + h2rX2 + 
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where: V(2, j) G {l,2,r}^, hij represents the channel gain of the link between nodes i and j, with the 
convention hrr = 0; Zi, and Z^. are additive white complex Gaussian noises with variances A''i, N2 
and A^r, respectively; the transmitted signals Xi, X2, X,. are subject to power constraints ]E|Xip < Pi, 
< P2 and E|XrP < -Pr- Note that, in Eq. ([T]), we dot not use any time index. This choice is made 
not only for the sake of clarity (time index is in fact required for proving the provided theorems) but also 
to indicate that we assume the propagation delays to be negligible. Always for the sake of clarity and 
simplicity, the channel gains between the different nodes (see Fig. [I]) will be considered fixed for the whole 
duration of the transmission. Physically speaking, this means that the assumed model can at least account 
for the path loss effects. In fact, by using standard arguments (see e.g., flU and ifTOl ). the presented results 
can be readily re-used for other types of channels like fast fading channels. Note that the expression of the 
signal transmitted by the relay (i.e., X^) depends on the protocol assumed and will therefore be provided 
in each of the three sections corresponding to the three types of protocols considered. What is general, 
however, is that the relay is assumed to operate in the full-duplex mode. Here again, this assumption has 
been well discussed in the related literature lfTT1l[fT2ll and will not prevent us from deriving results that are 
also applicable to the half-duplex mode, just as those obtained in other papers like lfT3l . 

Specific notations. Throughout this two-part paper we will use the following notations. The capacity 
function for complex signals is denoted by C{x) = log2(l + x). For any real a E [0, 1], the quantity a 
will stand for a = 1 — a. The notation —i will mean that — z = 1 if 2 = 2 and — i = 2 if z = 1. For any 
complex number, c G C, c*, |c| and Re(c) will stand for the complex conjugate, absolute value and real 
part, respectively. 

in. Transmission rates for the ZDSAF protocol 

In this section, the relay is assumed to implement an analog amplifier which does not introduce any delay 
on the relayed signal. The signal transmitted by the relay merely writes as Xr = a^l^ where corresponds 
to the relay amplification factor/gain. We will always assume a^. to be fixed. In the described setup, the 
following theorem provides a region of transmission rates that can be achieved when the transmitters send 
private messages to their respective receivers, the relay implements the ZDSAF protocol and the receivers 
implement single-user decoding (i.e., no multiuser detection or successive interference cancellation schemes 
are allowed here). 

Theorem 3.1 (Transmission rate region for the IRC with ZDSAF): Let Ri, i G {1,2}, be the Shannon 
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transmission rate for the source node Si. When ZDSAF is assumed the following region is achievable: 



where Pi = ^ and j = —i. 

The proof of this theorem is standard lfT4l and will therefore be omitted. The main point to be mentioned 
is that Gaussian codebooks have to be assumed to obtain the proposed region. An interesting point to 
investigate is the choice of the value of the amplification gain a^. In the vast majority of the papers available 
in the literature, is chosen in order to saturate the power constraint at the relay (E|Xrp = Pr) that is: 



can turn out to be sub-optimal in the sense of certain performance criteria. The intuitive reason for this 
is that the AF protocol not only amplifies the useful signal but also the received noise. While this effect 
can be negligible in certain scenarios for the standard relay channel, it is generally a dominant effect for 
the IRC. Indeed, even if the noise at the relay is negligible, the interference term for user i (i.e., the term 
hjrXj, j = —i) is generally not. This gives us a particular motivation for choosing the amplification factor 
ttr adequately that is, to maximize the transmission rate of a given user or the network sum-rate. The 
proposed derivation differs from lUSllfTTl because here we consider a different system (an IRC instead of a 
relay channel with no direct link), a specific relaying function (linear relaying functions instead of arbitrary 
functions) and a different performance metric (individual transmission rate and sum-rate instead of raw 
bit error rate [T5] and mutual information [[TtII ). Our problem is also different from ifTSl since we do not 
consider the optimal clipping threshold in the sense of the end-to-end distortion for frequency division relay 
channels. At last, the main difference with lfT6l is that, for the relay channel, the authors discuss the choice 
of the optimal amplification gain in terms of transmission rate for a vector AF protocol having a delay of 
at least one symbol duration; here we focus on a scalar AF protocol with no delay and a different system 
namely the IRC. In this setup, we have found an analytical expression for the best in the sense of Ri{ar) 
for a given user i E {1,2} and seen that the a^. maximizing the network sum-rate has to be computed 
numerically in general. The corresponding analytical result is stated in the following theorem. 

Theorem 3.2: [Optimal amplification gain for the ZDSAF in the IRC] The transmission rate of user i, 
Ri{ar), as a function of ar E [0,ar] can have several critical points which are the real solutions, denoted 



Vz E {1,2}, Ri<C 




(2) 




mentioned in some works [[T5llfT6ll[[T7l[fT8l . this choice 
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by a^j}} and afj, to the following equation: 

al [\mi\'^Re{piq*) - (jpip + Si)Re{min*)] + [|mi|2(|gi|2 + 1) - \ni\'^ {\pi\'^ + Si)] 

+ + l)Re(miO - |ni|2Re(pig*) = 

w/zere = hirKi^/pl, rii = ha^/pl, Pi = hj^Kiy/p], qi = hji^, Si = \hri\'^, i G {1,2} and j = -i. 

Thus, depending on the channel parameters, the optimal amplification gain a* = arg max Ri{ar) takes 

a value in the set a* G {0, a,., a^^j\ a^^^^}. If, additionally, the channel gains are reals then the two critical 
points write as: ai^ = and aj^ = _Z^i^?±!!™!^. 

The proof of this resuh is provided in Appendix |Al Of course, in practice, if the receive SINR (viewed 
from a given user) at the relay is low, choosing the amplification factor adequately does not solve the 
problem. It is well known that in real systems, a more efficient way to combat noise is to implement error 
correcting codes. This is one of the reasons why DF is also an important relaying protocol, especially for 
digital relay transceivers for which AF cannot be implemented in its standard form (see e.g., ifTSl for more 
details). 

IV. Transmission rates for the DF protocol 

The purpose of this section is essentially twofold. First, we state a corollary from [[U. Indeed, the given 
result corresponds to the special case of the rate region derived in where each source sends to its 
respective destination a private message only (and not both public and private messages as in [[T]|). The 
corresponding rate region is provided for making this paper sufficiently self-containing and for being used 
in the simulation part to establish a comparison between the different relaying protocols under consideration 
in this paper. Second, we explain why the DF protocol naturally introduces a game between the transmitters. 
The proof of the existence of Nash equilibrium in the corresponding game will be provided in Part II. The 
principle of the DF protocol is detailed in lfT3l and here we just give the main idea behind it. Consider 
a Gaussian relay channel for which the source-relay link has a better quality than the source-destination 
link. From each message intended for the destination, the source builds a coarse and a fine message. With 
these two messages, the source superposes two codewords. The rates associated with these codewords (or 
messages) are such that the relay can decode both of them reliably while the destination can only decode 
the coarse message. After decoding this message, the destination can subtract the corresponding signal and 
try to decode the fine message. To help the destination to do so, the relay cooperates with the source by 
sending some information about the fine message. Mathematically, this explanation translates in the IRC as 
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follows. The signal transmitted by Si is structured as = XiQ + ^J^^Xri where: the signals Xjo and 
X^i are independent and precisely correspond to the coarse and fine messages respectively; the parameter 
i>i represents the fraction of transmit power the relay allocates to user i, hence we have ui + V2 < 1; the 
parameter Tj represents the fraction of transmit power Si allocates to the cooperation signal. Using these 
notations we have the following result. 

Corollary 4.1 f/Ej/j.- When DF is assumed, the following region is achievable; for i G {1,2}, 

h,rWl-Ti)Pi \ ^ f \hii\''Pi + \Ki\''viPr + 2ne{huK^)^/nPiViPr \ \ 

Nr ; ' \ |%PP, + \hri\^ujPr + 27^e(%/i^) v^T.FjZ^jF, + j 
where j = —i, (z/i, z/2) G [0, 1]^ s.t. Ui + U2 < I and (ti, T2) G [0, 1]^. 

In a context of decentralized networks or in the case of unlicensed bands, each source Si has to optimize 
the parameter Tj in order to maximize its transmission rate Ri. As shown in the above rate region, this 
choice is not independent of the choice of the other source. Therefore, each source finds its optimal strategy 
by optimizing its rate w.r.t. r*(rj). In order to do that, each source has to make some assumptions on the 
value Tj used by the other source. This is precisely a non-cooperative game where each player makes some 
assumptions on the other player's behavior and maximizes its own utility. We therefore see that the DF 
protocol, through the parameter representing the cooperation degree between the source Si and the relay, 
introduces a PA game. In fact, one can even derive a hierarchical game since the relay can also be thought 
of as a player tuning (z/i, 1/2) (i.e., the fractions of power dedicated to users 1 and 2). All these issues will 
be treated properly in Part II. 

V. Transmission rates for the EF protocol 

In this section, we consider a third main class of relaying protocols, namely the estimate-and-forward 
protocol. The EF protocol can be implemented in digital transceivers (contrarily to the standard AF protocol) 
and always allows the receiver(s) to improve its (their) performance with respect to the non-cooperative 
case (in contrast with DF protocols that can degrade the performance in practice, see e.g., ifTSl ). In order to 
derive the corresponding transmission rate region for the Gaussian IRC, we have to prove that this region 
is effectively achievable. For this purpose, we first tackle the general case of discrete input discrete output 
channels. We know, from standard quantization and continuity arguments [[T4l . that the Gaussian case can 
be readily obtained from the discrete case. The proofs based on Shannon codes are not necessary for a 
good understanding and interpretation of the derived results and are therefore detailed in Appendices |B] and 
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O Before providing these results, we will first review the principle of the EF protocol which allows the 
reader to better understand the structure of this section. The principle of the EF protocol in the standard 
relay channel is that the relay sends an approximated version of its observation signal to the receiver. More 
precisely, in its information-theoretic standard version |fT3l , it consists in compressing the received signal at 
the relay in the Wyner-Ziv manner |fT9l i.e., knowing that the destination also receives a direct signal from 
the source that is correlated with the signal to be compressed. The compression rate is precisely tuned by 
taking into account this correlation degree and the quality of the relay-destination link. In our setup i.e., 
an IRC with two receivers, the EF protocol can be designed in at least two manners. On the one hand, the 
relay can define two resolution levels for its observation signal, which is what we call bi-level compression. 
On the other hand, it can also use a single resolution level adapted to the worse destination but reliably 
decodable by the better one (we call this scheme single-level compression). One of our contributions is to 
derive the transmission rates that can be achieved by using these two schemes and then discuss the difference 
between them. 

A. The case of discrete IRCs 

Before providing the two theorems associated with the two EF protocols proposed, we first define the 
discrete IRC. 

Definition 5.1: A two-user discrete memoryless interference relay channel (DMIRC) without feedback 
consists of three input alphabets Xi, and X^, and three output alphabets y^, 3^2 cind y^, and a probability 

n 

transition function that satisfies p{yi, 1/2, i/r \ x",X2, x") = JJ^p (?/i,fc, ^2,*:, ?/r,fc | xi^k,X2,k,Xr,k) for some 
n e N*. 

For the EF protocol based on a bi-level compression, the relay constructs two estimated versions of its ob- 
served signal Yr, which are denoted by Y^i and Yr2- As shown in AppendixO these estimates are constructed 
on a block-by-block basis: the relay constructs, from a block of n observations = (^/^(l), Urin)) two 
codewords and intended for Vi and V2 respectively. The signal transmitted by the relay results 
from mixing these two codewords. To this end, the coding scheme we used is superposition coding, which 
explains the presence of the coding auxiliary variables Ui, U2 in the theorem below; the meaning of these 
variables will be made obvious in the Gaussian case. 

Theorem 5.2: For the DMIRC {Xi x X2 x Xr,p {yi,y2,yr\xi, X2, Xr) , x 3^2 x X) with private mes- 
sages and bi-level compression EF protocol, any transmission rate pair (_Ri,i?2) is achievable, where 
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i?i < /(Xi; Fi, Yri\Ui) and R2 < /(X2; Y^, ^2^2), (5) 

under the constraints 

I{Yr;Yri\U,,Y,) < I{U,;Y,) and /(F,; ^21^/2, >^2) < I{U2;Y2), (6) 
/or some joint distribution 

P (Xi, X2, Ui, U2, Xr, Vl, ?/2, Vr, Vrl, yr2, Vrl, ?/r2) = 

p{Xi)p{x2)p (Ui) p {U2) P {Xr\Ui, U2) p (2/1, 2/2, VAxi, X2, Xr) P {yrllVr, Ul) P (?/r2|2/r, ^2) • 

In the case of single-level compression, the relay compresses Yr and encodes the same estimate Yr for both 
destination nodes. A natural question to be asked is whether the single-level resolution case is a special 
case of the bi-level case. The answer is no. Although the relaying principle is the same in both cases, 
as the Appendices show, there is one technical difference in the encoding procedure. In the bi-level case, 
the codewords used for superposition coding (associated with Ui and U2) at the relay are independent by 
construction. In the single-level case, there is a unique codeword for both destination. This is why making 
the choice Ui = U2 = in the bi-level case to obtain the single-level case is not admissible since the 
corresponding codewords would be totally correlated whereas they have to be independent. Technically, 
the choice Ui = U2 = Xr would violate the Markov chains assumed in Theorem 15. 2[ Rather, the case of 
single-level compression can be stated through the following theorem. 

Theorem 5.3: For the DMIRC {Xi x X2 x Xr,p {yi,y2,yr\xi, X2, Xr) , x 3^2 x X) with private mes- 
sages and single-level compression EF protocol, any transmission rate pair {Ri, R2) is achievable, where 

R,<l[x,-Y,X I ^r), ^e{l,2}, (7) 

under the constraint 

m.a.^I{Yr]Yr\Xr,Yi) < mmI{Xr]Yi), (8) 

i i 

for some joint distribution p{xi,X2, Xr, 2/2, 2/r, 2/ri, ?/r2, yr) = p(xi)p(x2)p(x^)p 2/2, yA^i, X2, Xr) p (yr|?/r, Xr) . 

B. The case of Gaussian IRCs 

As mentioned in the beginning of this section, obtaining achievable transmission rates for Gaussian 
IRCs from those for discrete IRCs is an easy task. Indeed, the latter consists in using Gaussian codebooks 
everywhere and choosing the coding auxiliary variables properly (i.e., Ui, U2, Yr^i and Yr,2 in the bi-level 
case and, Xr and Yr in the single-level case). 
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For the bi-level compression, the coding auxiliary variables Ui and U2 are chosen to be independent and 
distributed as Ui ~ A/'(0, uiPr) and U2 ~ A/'(0, 2^2 ^r)- The corresponding codewords and convey the 
messages resulting from the compression of Yr. The auxiliary variables F^. 1, Yr,2 write as Yr,i = Y^ + Zw} and 
Yr,2 = Yr + Zwl where the compression noises Zwl ~ Ar(0,A^4'i) and Zw\ ~ N'{0, Nw}) are independent. 
At last, the relay transmits the signal = U1 + U2 as in the case of a broadcast channel except that, here, 
each destination also receives two direct signals from the source nodes. By making these choices of random 
variables in Theorem 15.21 we obtain the following theorem. 

Theorem 5.4: For the Gaussian IRC with private messages and bi-level compression EF protocol, any 
rate pair {Ri, R2) is achievable where 



1) if c 



|ferl|'i^2Pr 



|/lllPPl + 1/121 PPa + \hrl\^ViPr + N-L 
( 

Ri < C 



> c 



\hr2\ 



\hn\'Pi 



\h22\^P2 + \hl2\^Pl + \hr2\^UiPr + N2 



, we have 



+ 



I 



Ri < C 



\h2r-\'^P2+Nr- + Ni^} 

\h22?P2 



\h2r?P2 



\ 



\N2 + \hr2\''yiPr 



\hlr\^Pl+Nr+N, 



(9) 



(10) 



subject to the constraints N^} > ('^"'^^^'^i;^;!;!^^-^^-^^ and N^} > 



2) else, if C 



\hr2\ VlPr 



|/122PP2 + |/ll2|2Pl + \hr2\'V2Pr + N2 
( 

Ri < C 



hrll^VlPy 

> c 



I 1 I ^1 -Pr 



\h^2\'^U2Pr 



\hn\'Pi 



1/illpPl + |/l21pPl + \hrl\^U2Pr + Nl 
\hlr\^Pl 



, we have 
\ 



\ho^\^P2(N^ + Ni}}) a; , AtW , \h2r\^P2{\h^l\^l^2Pr+Nl) 

\ Nl + \hrl\^V2Pr + ^ ^ + ^^'"^ + \h2i\-P2 + \h.l\-^2Pr+m 



R2 < c 



\h22?P2 



\h2r?P2 



i/..2pp.K+A-i-i) Nr + N^} + i:^;;^^ 

\hlr\^Pl+Nr+N^^ 



(11) 



(12) 



subject to the constraints > > (|/^..|-p.+I^..|-Pi+^.)^-^^ 



l/iril^i'iP. 



|/lrfpl'2Pr 



3) 



Ri < C 



R2 < C 



\hii\^Pi 



\hlr\^Pl 



\ Nl + \hrl\'^V2Pr . 

\ ^ ^ I '^-^l Ih2r.\^p^+N^ + N.. 

\h22?P2 



'ii2PPifA^r+JVil'') /V J- /v(^) -I- 



\hlA^Pl{\hr2\'^VlPr- + N2) 
\hl2\'^Pl + \h^2\'^VlPT+N2 



(13) 



(14) 
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subject to the constraints nL'} > i t^l^/X"'^^^''^^^'^^ ^ (|...pp.+|.,.|^P,+|.^^VP.+Ar.)^-^i 

with {iyi,iy2) e [0, + 1^2 < I, A = |/ii,.pPi + |/i2rP-f2 + Nr, Ai = 2'Re{hiihl^)Pi + 27^e(/l2l/l;JP2 ^2 = 

27^e(/^l2/^^ jPi + 27^e(;^22/^2r)^'2■ 

The three scenarios emphasized in this theorem corresponds to the following three situations: \)Vi has the 
better link (in the sense of the theorem) and can decode both the relay message intended for V2 and its 
own message and can therefore cancel the interference signal due to the relay; 2) this scenario is the dual of 
scenario 1); 3) in this latter scenario, none of the destination is able to suppress the interference generated 
by the relay. 

For the single-level compression, one just needs to choose YJ, as 1^ = 1^ + with ~ A/'(0, A'^^) 
and apply Theorem 15.31 We find the following result. 

Theorem 5.5: For the Gaussian IRC with private messages and single-level compression EF protocol, the 
rate pair (i?i,/?2) is achievable, where 

T) I I ^^^^ I "^"^ I ^^iV I I / -t ^\ 

P-i — ^ \ ~ \h ..\2p.tM :M T "I" ~~ ~~ \h . 12 p. AT. (15) 



\hj,\^Pj{Nr+N^.) ' AT , AT , IfejrPPjjV, 



max-! tT?,(T2 f 

subject to the constraint N^z > — ^2^0 , with 



Ro = min C „^ '77 , ,C ,, '77 ; „ , , (16) 





\hrl\ 


1 J. 


\hu\ 


?Pi + 


\h2l\ 


PP2 + iVl 





\hr2\ 


\2p 




\'P2 


+ 


\hi2\ 


PPl+iV2 



-I. 



J 2 II, I2n , II, I2n , Ar (2Tle(h,ihl)Pi+2-Re(h,,h* )PjY ■ ^ n oT j • 

With this type of EF protocols, the relay adapts the compression resolution to the worse receiver. At this 
point, one interesting question arises. Is it better to choose the bi-level compression EF protocol or the single- 
level compression EF protocol? This is one of the issues treated in the next section, which is dedicated to 
simulations. It turns out that, in some useful special cases, an analytical answer can be given. This is the 
purpose of the end of this section. 

The first special case is the case of asymmetric channels. By asymmetric we mean that ^ A'^i or 
A^i ^ N2. Consider w.l.o.g. the case N2 ^ A^i with N2 +00 and A^^i < +00. It is obvious and easy to 
check that R2 ^ for both the single-level (SL) and bi-level (BL) compression schemes. For the rate of 
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user 1 we have: 

?(SL) 





hn 


'Pi 


\h2i 


VP2 + A^i 



lim B}r' = C , 'V, ' ^ (17) 



lim i?S^^^ = C 

N2^oo 



( \ 

J- ' 17. lOn iTVT lAr(l) 



(18) 



It is clear that user 1 gets a higher rate with the bi-level compression scheme and, as the rate for user 2 
tends to zero, this protocol is also the better one in terms of network sum-rate. A second case, which is 
very special but has the advantage of requiring no long derivations is the perfectly symmetric IRC, which 
we define by Pi = P2, Ni = N2, \hir\ = \h2r\, \hi2\ = \h2i\, |/iii| = |/i22| and \hri\ = \K2\- With this 
setting, a quick inspection of the compression noise expressions in Theorems 15.21 and 15.31 shows that the 
compression noise level (i.e., A^„,^) of the single-level scheme will always be lower than those (i.e., N^I 
and Nw}) obtained with the bi-level scheme. This is because, in symmetric IRCs, in addition to the fact that 
both receivers can perform almost similarly in terms of wyner-ziv compression noise level, the single-level 
compression scheme does not generate additional interference whereas, in the case of bi-level compression, 
there is always at least one user undergoing some additional interference from the relay. 

VI. Simulation results 

Simulation setup. In all this section we assess the performance of the IRC in terms of system sum-rate as 
a function of the relay position. For this purpose we assume a path loss model for the channel gains \hij\ 
and a given location for each node. For the path loss model we take \hij\ = (^^^ ^ for E {l,2,r}^ 
where do = 5 m is a reference distance and 7 = 2 is the path loss exponent. The nodes Si, S2, Vi, V2 are 
assumed to be in a plane. The positions of the nodes will be indicated on each figure considered and are 
characterized in this plane by the distance between the nodes which are chosen as follows: d'n = 11.5 m, 
£^22 = 10 m, (i'i2 = 11 m and d'21 = 14 m. As for the relay, to avoid any divergence for the path loss in 
dij = 0, we assume that it is not in this plane but at e = 0.1 m from it i.e., the relay location is given by 



the {xr, Vr, Zr) whcrc Zj. is fixed and equals 0.1 m; thus d^ = y d'^j + for i = r or j = r and i j. The 
noise levels at the receiver nodes are assumed to be normalized and unitary (iVi = N2 = Nr = 1). In terms 
of transmit power we analyze two cases: a symmetric case where Pi = P2 = 10 (normalized power) and 
an asymmetric one where Pi = 3 and P2 = 10. The relay transmit power is fixed: = 10. We have seen 
that, for both DF and bi-level compression EF protocols, the relay has to allocate its power between the two 
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cooperative signals intended for the two receivers. We consider two different cases: the uniform PA policy 

(vi = P2 = \) and the optimal PA policy (z/J', vl) = argmax(,^j^j,2) Ri + R2 s.t. z^i + 1^2 < 1 that maximizes 
the system sum-rate. 

Comparing the dijferent relaying protocols. We start by investigating the best choice for the relaying 
protocol between the AF, DF and bi-level compression EF protocols. For the symmetric scenario, Fig. |2] 
represents the regions of the plane G [—4, +4] x [—3, +4] (corresponding to the possible relay 

positions) where one given protocol performs better than the two others in terms of system sum-rate. These 
regions are in agreement with what is generally observed for the standard relay channel. This type of 
information is useful, for example, when the relay has to be in some places because of some practical 
constraints and one has to choose the best protocol for a given location. Fig. |3] allows one to better quantify 
the differences in terms of sum-rate between the AF, DF and bi-level EF protocols since it represents the 
sum-rate versus Xr for a given j/r = O.Sdo- The discontinuity observed stems from the fact that for the 
bi-level EF protocol there is a frontier delineating the scenarios where one receiver is better than the other 
and can therefore suppress the interference of the relay (as explained in Sec. IV-B|) . We indicate this frontier 
in the figures commented just below. 

Single level compression vs. bi-level compression. Here we focus on the EF protocol and its two variants 
proposed. Fig. |4] corresponds to the symmetric scenario and = 1^2 = \ and shows the regions where 
the relay performs better by using the bi-level compression scheme instead of the single-level one. The 
dotted line represents, for the bi-level compression EF protocol, the frontier separating the two scenarios 
"receiver 1 is better than receiver 2" and conversely (in the sense of Theorem 15.41) . These regions illustrate 
a combination of several trade-offs, which we try to explain clearly here. When the qualities of the links 
IZ ^ Vi and IZ Vj (with j = —i) are quite similar and it is also the case between the links Si ^ IZ 
and Sj TZ, the single-level scheme is the most efficient because: the compression rate for each receiver 
is not limited by the other receiver; the relay uses all its power Pj. to transmit the cooperative signal; as 
there is only one cooperation message, the relay generates no additional interference. This scheme is thus 
the best choice both for Ri and R2 and therefore for Ri + R2. On the other hand, if the qualities of the links 
IZ ^ Vi and IZ Vj are markedly different, which happens when the relay is located around Vi (resp. 
Vj), it is better for Ri (resp. Rj) to choose the double-level scheme. This is true if the negative effect that 
the double-level scheme allocates only -y to each cooperation message is compensated by the positive effect 
that the compression noise level seen at Vi (resp. Vj) will no longer be lower-bounded by the one at Vj 
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(resp. Vi). Fig. |4] precisely illustrates the complex trade-offs just mentioned. Fig. |5] shows these trade-offs 
for the asymmetric scenario. 



In this part we have derived achievable transmission rate regions for the interference relay channel when 
the relay respectively implements the AF, DF and EF protocols. In comparison with the standard relay 
channel, we have emphasized several important differences. For the AF protocol, from a given receiver's 
point of view (say from Vi), the relay not only amplifies the noise but also the interference due to ^2. 
As a consequence, scenarios where the amplification factor used by the relay does not saturate the power 
constraint at the relay are much more frequent than for a relay channel. We have seen that the optimal 
factor can be determined analytically. For the DF protocol we have seen to what extent it creates a game in 
interference relay channels. For the EF protocol, we have seen that because of presence of several receivers, 
the EF protocol can be designed based on multi-level compression, which leads to several variants of the 
EF protocol (the single-level and bi-level compression schemes in our case of two receivers). The best 
choice between them depends on the performance criterion to be optimized (a given individual rate or the 
network sum-rate) and the network topology (symmetric/asymmetric networks). The provided theoretical and 
simulation results gives insights on how to choose the relaying protocol in IRCs or locate a relay optimally 
for a given protocol. An extension of this work is proposed in Part II where we investigate the case of 
multi-band IRCs with cognitive transmitters that can freely and selfishly allocate their power between the 
available bands. 



Using the notations given in Theorem 13.21 and also the SINK in the capacity function of Eq. (|3.1|) the 
rate Ri can be written as: 



We observe that Ri{0) = C y \g"\i+i j that we have an horizontal asymptote 

(I Tn ■ P \ 
-. — 777^ I . Also the first derivative w.r.t. Ur is 
iPir + SiJ 

. _al [|m,|^Re(p,i7*) - {\pi\^ + Si)Re(mi?i*)] + [|m,|^(|qi|^ + 1) - k»P(|p«P + Si)] + + l)Re(min*) - |nipRe(p,g*) 

[\Piar + <?iP + Sial + l][|miar + 7i,j2 + Ip^a,. + g^P + Sjffl? + 1] 



VII. Conclusion 



Appendix A 



Optimal amplification gain at the relay with the ZDSAF protocol 
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The explicit solution, a* depends on the channel parameters and is given here below. We denote by 
A the discriminant of the nominator in the previous equation. If A < 0, then in function of the sign of 
\mi\'^Re{piq*) — + Sj)Re(mj?7,*), the function Ri{ar) is either decreasing and a* = or increasing and 
a* = a^. Let us now focus on the case where A > 0. 

1) If \mi\'^Re{piq*) - + Si)Re{min*) > then 

a) if a'-^j < and a'-^j < then a* = ar; 

b) if aJJ > and aJJ < then 

i) if ar > then a; = 0; 

ii) if Hr < o-i^i then 

• if -Ri(O) > Riittr) then a* = else a* = 

c) if a^j^} < and a^^} > then the analysis is similar to the previous case and a* G {0,0,} 

(2) 

depending on Or this time; 

d) if > and aJJ > 

•^ -f (1) ^ (2) 

i) if a^J < a).i 

A) if Hr < a'l^J then a* = a^; 

B) if aj.^' < ar < aj.5 then a^ = aJJ; 

C) if ttr > afj then 

. if Ri{a^^}) > Ri{ar) then a* = a|.^/ else a* = ar', 

ii) if a^^j > af] then the analysis is similar to the previous case, exchanging the roles of a!-^} 
and Or^^; 

iii) if a'^J = a^^J then a* = ar. 

2) If \mi\'^Re{piq*) - + Si)Re{min*) < then 

a) if aJ5 < and aJJ < then a; = 0; 

b) if aj.^^ > and aj.^ < then 

i) if ttr > air,} then a* = a,, else a* = a^^} ; 

c) if a^j^j < and a^j^j > then the analysis is similar to the previous case and a* G {ar,al^-} 
depending on a^j^j this time; 

d) if > and aj.} > 

•\ -f (1) ^ (2) 

1) if alj < alj 
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A) if ttr < a^^i then = 0; 

B) if a^^i <ar< aij then 

. if Ri{0) > Ri{ar) then a* = else a* = Ur, 

C) if ttr > afj then 

. if Ri{af}) > Ri{0) then a*, = 4J else a; = 0; 

ii) if a^j^j > afj then the analysis is similar to the previous case, exchanging the roles of a:^J 
and af^; 

iii) if o'^J = afj then a* = 0. 

We observe that, depending on the channel parameters, it is not always optimal to saturate the power 
constraint at the relay. 

Appendix B 

achievability proof of the rate region in theorem [5^ 

Definitions and notations 

We denote by A^^'' {X) the weakly e-typical set for the random variable X. If X is a discrete variable, 
X ^ X, then denotes the cardinality of the finite set X. We use to indicate the vector (xi, 2:2, ... , 

Definition B.l: A (2"^\ 2"^^ , n) -co for the DMIRC with private messages consists of two sets of 
integers Wi = {l, 2"-^^! } and W2 = {l, ...,2''^^}, two encoders: fi : Wi X^",, a set of relay 
functions {fr,k}k=i ^^^^ ^^^^ ^r,k = fr,k {yr,i,yr,2, ■■■,yr,k-i) , 1 < k < 71 and two decoding functions 
di '■ y'i' ^ « G {1,2}. The source node Si intends to transmit Wi, the private message, to the receiver 
node Vi. 

Definition B.l: The average probability of error is defined as the probability that the decoded message 
pair differs from the transmitted message pair; that is, Pi"^ = Pr [gi (F") 7^ Wi or g2 (F2") 7^ ^2 | (W^i, W2) ] , 
where (^1,^2) i^ assumed to be uniformly distributed over Wi x W2. We also define the the average 
probability of error for each receiver as Pl:f^ = Pr [gi (YJ") 7^ Wi \ Wi] . We have < max |-Pi"\ Pe2^^ — 
Pe^'' < Pg"^ + Pg^''. Hence Pi"^ implies that both P^f ■* and P^2^ 0, and conversely. 

Definition B.3: A rate pair (^1,^2) is said to be achievable for the IRC if there exists a sequence of 
(2"-f^i, 2"^^ n) codes with Pi"^ ^ Q as n ^ 00. 

Overview of coding strategy 
At the end of the block k, the relay constructs two estimations y"i(A;) and y"2(^) of its observation 
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that intends to transmit to the receivers Vi and 1^2 to help them resolve the uncertainty on wi^k and W2,k 
respectively at the end of the block A; + 1. 
Details of the coding strategy 

Codebook generation 

i Generate 2'"^' i.i.d. codewords x^{wi) ~ 11^=1 

where Wi e {l, . . . , 2"-^"}, i e {1,2}. 

ii Generate 2"'^o ^ i.i.d. codewords w" ~ 11^=1^ ("^i.fc)- Label these Si G |l, . . . , 2"'^o ^ 1. 

iii Generate 2"^o i.i.d. codewords U2 ~ YYk=i P i'^2,k) ■ Label these ^2(52), si e < 1, . . . , 2"^o >. 

iv For each pair (m"(si), (S2)), choose a sequence where x" ~ p (S2)) = 

Yll=lP{Xr,k\^l,k{^l)^ U2,fc(s2))- 

V For each ^"(si), generate 2"^^ conditionally i.i.d. codewords y"]^ ~ 11^=1 (^/f ifcki.fcl'Si)) and 
label them yl^^{zi\si), zi E 1 1, . . . , 2"^^ |. For each pair {ui,yri) E Ui x y^-i, the conditional 

probability p{yri\ui) is defined as p{yri\ui) = ^ p{xi)p{x2)p{yi,y2,yr\xi,x2,Xr)p{yri\yr,ui). 

x-i,x-2,yi,y-2,yT 

vi For each ^2(52), generate 2"^^ conditionally i.i.d. codewords ~ 11^=1 P {yr2k\u2,k{s2)) and label 
them y"2(^2|s2), 2^2 G |l,2"'^2| each triplet {u2,yri) EU2 x yri, the conditional probability 

Piyr2\u2) is defined as p(2/r2|M2) = ^ p (Xi) p (X2) P ?/2, ?/r|Xl, X2) P (2/r2|Z/r, ^2)- 



xi,x2,yi,y2,yr 



vii Randomly partition the message set |l, 2, . . . , 2""^^^ into 2"^o'^ sets js'f \ 5*2^^ . . . , | by 
independently and uniformly assigning each message in 1 1, ... , 2"^^ | to an index in 1 1, ... , 2"-^o ^ | 

viii Also, randomly partition the message set 1 1 , 2, . . . , 2"^^'' | into 2'^^o' sets | Sj^^ , S'f \ . . . , | 
by independently and uniformly assigning each message in 1 1 , . . . , 2"^^^ | to an index in 1 1 , . . . , 2"^o 

Encoding procedure Let wi^k and W2,k be the messages to be send on block k. Si and ^2 respectively 
transmit the codewords and X2(w2,fc)- We assume that {ui{si^k-i) , yriizi,k-i\si,k-i) , yrik — 1)) E 

A'f'^ and zi^k-i ^ •S'ij-'j^ and also that 

{U^{s2,k-l),yr2iz2,k-l\s2,k-l),yrik - 1)) E with 

Z2,k-i G 'S'i2,l- Then the relay transmits the codeword 

<(Sl,fc,S2,fc)- 

Decoding procedure In what follows, we will only detail the decoding procedure at the receiver node Vi 
(at V2 the decoding is analogous). At the end of block k: 

i The receiver node Vi estimates Si ^ = Si if and only if there exists a unique sequence that 
is jointly typical with yi{k). We have si = si ^ with arbitrarily low probability of error if n is 
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sufficiently large and R^(^^ < I{Ui;Yi). 

ii Next, the receiver node Vi constructs a set Ci (y^ {k — 1)) of indexes zi such that 

{ui (si^k-i) , Vri (^il-si.fc-i) ,yi {k — 1)) G Ai^K Vi estimates zi^k-i by doing the intersection of 
sets Ci{yi{k — 1)) and 5"^^^^. Similarly to |fT3] theorem 6] and using |fT3l, lemma 3], one can 
show that 5i = -2i,fc-i with arbitrarily low probability of error if n is sufficiently large and 
Ri<I{Yri;Yi\Ui) + R^^\ 

iii Using and y"(A; — 1), the receiver node Vi finally estimates the message wi^k-i = 
Wi if and only if there exists a unique codeword x^{wi) such that 

{x'l{wi),u'^{si,k-i),yi{i - l),yri(-Si,fc-i|si,fc-i)) e A'e'\ Ouc cau show that wi = wi^k-i with 
arbitrarily low probability of error if n is sufficiently large and 

Ri < I (^Xi;Yi,Yri I f/i). (19) 

iv At the end of the block k, the relay looks for the suitable estimation of its observation that it 
intends to transmit to the receiver node Vi by estimating fi fc. It estimates 5i ^ = zi if there exists 
a sequence such that {ui{si^k) , yri{zi\si^k) , yr{k)) G Ai^\ There exists a such sequence 
if n is sufficiently large and Ri > I{Yri; Yr\Ui). 

From i, ii, iii we further obtain 

I{Yri;Yr\Ui,Yi) < I{Ui;Yi). (20) 

(fT9l) and (l20l) are identical to the first terms in ([5]) and respectively. The achievability proof for the 
second receiver node follows in a similar manner. Therefore, we have completed the proof. 

Appendix C 

Achievability proof of the rate region in Theorem [53] 

Details of the coding strategy 

Codebook generation 

i Generate 2'"^' i.i.d. codewords x'^{wi) ~ Y[k=iP where Wi e {l, 2, . . . , 2"-^'j, i E {1, 2}. 

ii Generate 2"^" i.i.d. codewords xj! ~ YYk=iPi^r,k) and label them xj!(s), s G {1,2"-^°}. 

iii For each codeword x"(s), generate 2"^ conditionally independent codewords ~ 11^=1 P iyr,k\xr,k{s 
and label them y"(z|s), z G {l,...,2"^}. For each pair {xr,yr) G x y^, the conditional 

probability p{yr\xr) is defined as p(y.r|.x,.) = ^ p{xi)pix2)p{yi,y2,yr\xi,x2,Xr)p{yr\yr,Xr)- 
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iv Randomly partition the message set { 1, 2, . . . , 2"^} into 2"^" sets {5*1, 5*2, ... , S2nRo } by indepen- 
dently and uniformly assigning each message integer in |l, . . . , 2"^} to an index in |l, . . . , 2"^"}. 

Encoding procedure Let wi^k and W2,k be the messages to send in block k. Si and S2 respectively 
transmit the codewords x'^{wi^k) and X2{w2,k)- We assume that {x^{sk^i),y^{zk^i\sk^i),y'^{k — 1)) G Ai^^ 
and Zk-i G Ss^^. Then the relay transmits the codeword xl'{sk). 

Decoding procedure The decoding procedures at the end of block k are as it follows: 

i As a first step, the receiver node Di, i E {1,2}, estimates = s if and only if there exist a 
unique sequence x^{k) that is jointly typical with y^{k). We have s = Sk with arbitrarily low 
probability of error if n is sufficiently large and Rq < /(X^;!^), Vi G {1,2}, which we rewrite 
as Ro < min{/ (X,; Y^) , / (X,; Y2)}. 

ii Next, the receiver node constructs a set £j (y" {k — 1)) of indexes z such that 

(x" (sfc_i) , (2;|sfe_i) , (A; — 1)) G Ai"\ Vi estimate by doing the intersection of sets 
Ci (y" {k — 1)) and Sg^^. Similarly to [[III theorem 6] and using [[131 lemma 3], one can show that 
Zk-i = Zk-i with arbitrarily low probability of error if n is sufficiently large and R < KYr'-, Yi\Xr) + 
Ro, Mi G {1, 2}, which can be written as i? < min |/ {Vr] Yi\X^ , / [Vr] Y2\X.^ j + Rq. 

iii Using y'^{zk_i\sk-i) and ?/"(A; — 1), 2 G {1, 2}, the receiver node Vi finally estimates the message 
Wi^k-i = Wi if and only if there exists a unique codeword such that 

(xf (wi), 2/r(A;-l), y"(4-i|Sfc-i)) G At\ One can show that Wi = Wi^k-i with arbitrarily 

low probability of error if n is sufficiently large and 

R^ < l[Xf,Y„Yr I X,). (21) 

iv At the end of the block k, the relay looks for the suitable estimation of its observation by estimating 
Zk- It estimate Zk = z if there exists a sequence yl!-{z\sk) such that {x'^{sk), y^{z\sk), yr{k)) G Ai"\ 
There exists a such sequence if n is sufficiently large and R > I{Yr', YA^r)- 

From ii and iv we have that /(i;;y^|X^) < min |/ (Yr;Yi\Xr^ , / (^i;; K^l^r) }+min {/ (X^; Fi) , / (X^; Ys)} 
which is equivalent to 

max|/(y,;i;|X,,yi),/(i;;F,,|X„y2)} < min{/(X,;ri),/(X,;F2)}, (22) 

because of the Markov chain Y^ Q {Xr,Yr) Q Yi, Mi G {1,2}. ([IB and ^ are identical to © and ([8]), 
respectively. 
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Fig. 2. For different relay positions in the plane £ [—4, +4] x [—3, +4], the figure indicates the regions where one relaying protocol 

(AF, DF or bi-level EF) dominates the two others in terms of network sum-rate. 
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Fig. 3. Achievable system sum-rate versus x,- (abscissa for tlie relay position) for a fixed j/r (j/r = 0.5do), with AF, DF and bi-level EF. 
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X / d 



Fig. 4. Symmetric scenario in terms of transmit power: for different relay positions in the plane ^f^i j G [^2, +3] x [—2, +3], the figure 
indicates the regions (around Vi) where the bi-level compression EF protocol performs better than its single-level counterpart in terms of 
network sum-rate. 
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Fig. 5. Asymmetric scenario in terms of transmit power: for different relay positions in the plane (^f^, G [—2, +3] x [—2, +3], the 
figure indicates the regions (around Di) where the bi-level compression EF protocol performs better than its single-level counterpart in terms 
of network sum-rate. 
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